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Abstract : In a previous paper jTJ, we have shown that the gamma subordi- 
nators may be represented as inverse local times of certain diffusions. In the 

present paper, we give such representations for other subordinators whose 

C 

Levy densities are of the form - — , . , < 7 < 2, and the more general 
_ _ (smh(y))j' 

family obtained from those by exponential tilting. 



Keywords : Subordinators, Krein correspondence, inverse local times. 



1 



1 Aim of the paper and summary of [1] 



1.1 In this paper, we continue the program started in [T], that is : to 
represent as many subordinators (Se, £ > 0), i.e : increasing Levy processes, 
started at 0, as possible as inverse local times (ti, £ > 0) of some particular 
M + -valued diffusion (X t ), such that is regular for itself (relatively to X). 
More precisely, assume that : 

E[exp(-XS e )} = exp(-^(A)) , 

POO 

where : ^(A) = / v{dy){l — e~ Xy ) } and v{dy) - the Levy measure associated 
Jo 

with (iS>£, I > 0) - is of the form : 

POO 

v{dy) = h{y)dy , with h(y) = / d\(x) e~ xy , 

Jo 

for some positive a-finite measure \(dx) on R + , then, it is known, as a 
consequence of Krein's theory (cf : Knight [5], Kotani-Watanabe that 
there exists a unique diffusion (X t ) taking values in K + , such that its inverse 
local time at 0, (17, £ > 0), is distributed as (Se, £ > 0). 
Finding X when v is given is called (here) Krein representation problem. 
In our paper flj, we could fill in the following 

Table 1 



%) 


Generator of (X t ) 


Distribution 


c 

(0 < a < 1) 


^ = l£ + ^£; * = 2(l-a) 




ya +i e 
(0 < a < l;fi > 0) 


where ^(2/) = y a K a (y), y>0 


p<5;/4 


0*>o) 


& = § A + (± + v^g £ 


p2;/4 
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1.2 Details of Table 1 



In fact, the result for the first row goes back at least to Molchanov-Ostrovski 
|H], the result for the second row is deduced from that in the first row with 
the help of the following discussion, which relates Esscher transforms (of 
subordinators) to Girsanov transforms (of diffusions). 
If (X 4 ) t > is a M + -valued diffusion, whose inverse local time at : 

r e = inf{t : L t >£},£> , 

admits Levy measure v{dy), and Levy exponent (^f(d),9 > 0), and if one 
defines : 

<p gi {x)=E x [exp{-0T o (X))] 1 

then, there is another diffusion, which we shall denote by (X t \t > 0), with 
laws (P^,x > 0), such that : 

ni=^§-MmL t -et).p^ (i.i) 

whose inverse local time {r^l > 0) under P^ satisfies : 

E% l (exp -\T t ) = exp(-£(^(X + 6)- V(0))) 

i.e : this inverse local time is the ^-Esscher transform of (re, £ > 0) under 
P : its Levy measure (under P^) is : e~ 6y u(dy). 

It is also noteworthy that, under some adequate restriction of their domains, 
the infinitesimal generators C 9 ^ and £ are related by : 

£»> = £ + A(lo g(mW )).|- 

Finally, the result for the third row was deduced by letting a — > in the 
second row, while taking care of the choice 1 of the local times made for 
C- a> iJ,i- (A compendium of choices of local times for Bessel-like diffusions is 
made in |2J). 

1.3 In the present paper, we wish to complete the preceding Table 1, by 
considering the 3 parameter family of Levy measures on M + : 

( \ a+l 
sinh(/xy) ) exp ^) dy ( L2 ) 



x As is well-known, the local time in a standard Markovian set up, at a given level, is 
unique up to a multiplicative constant, which for our studies, needs to be chosen carefully. 
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(The "true" parameters are : fi > 0, k, and a; as before, C is simply there 

to ensure an additional degree of freedom, if necessary). 

In order that Vfx ja ,k{dy) be a Levy measure, i.e : it must satisfy 

POO 

/ (xAl) V/j,,ai,k (dx) < oo , we need : 0<a<l;A;<l + a. 
Jo 

We now recall that, from Pitman- Yor |TT| formulae (16), p. 276), if Q 5 Z ' M , for 
< 5 = 2(1 — a) < 2, and fi > 0, denotes the distribution of the squared 
radial Ornstein-Uhlenbeck process, with "dimension" 5, and parameter //, 
started from z, i.e : the solution of : 

dZ t = 2^Z~ t dB t + (5 - 2fiZ t )dt ; Z t >0,Z = z, 

then, under Qo^) the inverse local time (re,£ > 0) admits as its Levy 
measure : ^ 

C (^^)) eXp(/i ^ )rf2/ (L3) 

which is a particular case of (I1.2L with k = - = (1 — a). 

In the next section, we shall show, essentially with the help of the recipe 
1)1. lj) . how to construct a diffusion, indexed by the 3 parameters (a,/j,,k), 
which solves the Krein representation problem for v^ a ^. 

1.4 Some among the new diffusions we are finding as solutions of Krein's 
problem are related to the diffusions we found in |T] by time changing. We 
first discovered this relationship by applying the analytical identity : 

Wo,p(z) = ^Kp (|) (1.4) 

between Wo,/3 , a Whittaker function with parameters (0, (3), and Kp (see 
Appendix). Thus, a part of our present discussion may be considered as 
giving a probabilistic interpretation to 1)1.4)1 . 
We also develop a similar discussion for the analytical identity 

M ,^) = 4^r(/3 + l)V5^(|). 

2 Solving Krein's problem for v^ a ,k 

2.1 We take up the notation in ()1.3j) : in fact, it is more convenient to 
consider the family of radial Ornstein-Uhlenbeck processes (and not their 
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squares), which we shall denote as {R 5 ^{t), t > 0) and their laws (P^, 
r > 0). It will be helpful, for the sequel, to have the following formula at 
hand, for the infinitesimal generator C- a ,n of R 5,fl : 

Id 2 [5—1 \ d 

It is well-known (see, e.g., Pitman- Yor |9J, p. 454, formula (6.b)) that there 
is the relationship : 

/ 2(d _ i \ 

B*>f(t) = e-^R s I 2 J > t > , ( 2 - 2 ) 

where, on the RHS, (R s (u), u > 0) denotes a 5-dimensional Bessel process. 
Thus, we obtain : 

' e ^;P^) ( = l (r„;P I i ) < '= ) ^ (2.3) 

ZjJ J (a) (b) 2-f a 

where, on the RHS, 7 a denotes a gamma variable with parameter a. [(a) 
follows from (|2.2j) . while (b) is well-known, and goes back to Getoor PJ, see, 
e.g., Yor [13], for some variants...]. 

We thus deduce the following formula from (12. Hj) . with the help of elementary 
computations : 

JS5"(exp(-«T )) = — -r / " — - dt 



r(a) 



(/xx 2 ) — e"TW (1 _ a) _ A (/ix 2 ) for -2a/i<# (2.4) 
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where denotes the Whittaker function, with parameters (a, b). 
2.2 We now write : 

( \ Q+l 
sinh(/xy) ) e Mv k y) d y 

( \ Q + l j 
-j^yj exp( / ,-,)exp(-^)^, 

where : = // ^— — kj , k < 1 + a. 

According to the preceding computation, we now find that the diffusion with 
infinitesimal generator : 

= + ^ i g { (/ ^)^V^V (1 ^ ^ >x 2 )} .± (2.5) 

5 



solves Krein's representation problem for z/ M a (We note in fact that : 
(1 — a) - - k 

^ = -, so that : W (1 _ a) _ % JO = ^.#(0)- 
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The case where k = is particularly interesting, since, on one hand 

a+l 



v»,a,o{dy) = C 



IX 



sinh(/xy) 

and on the other hand (see Appendix) : 



dy, 



7T 2 V2 



(2.6) 



so that the diffusion which solves Krein's representation problem for v^a^dy) 
is the solution to : 



dX t = dB t 



6-1 
~2X7 



'K'« 



dt 



Here, we need to give some details about this computation : 

We deduce from formula (I2.5jl . in the 
with the help of formula 1)2.611 . that : 



(2.7) 



5ll 



a) We deduce from formula (I2.5jl . in the particular case fc = 0, i.e : # = — 



b) Now, trivially : 



and, equally simply 



a 



- "* *» (/*) 



En 



which translates into the stochastic differential equation form of for- 
mula (ED) . 
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As an introduction to the next discussion, we write down 

Table 2 



%) 


Generator of (X t ) 


Distribution 


ya +i e 
(0 < a < l;/! > 0) 




(5 = 2(1 -a)) 


n ( m \ a+1 ^s y 

ysinh(fiy) J ' 

0<<y = 2(l-ot) <2;//>0 






/ m \ a+1 
V sinh (^) / 

The first row is simply taken 


^-a,n — 1 — a T ftx \ V> 2 ) dx 

from Table 1 (second row there). 





As said above, the second row follows from Pitman- Yor [TT]. In the third 
row, we have written C_ a for the infinitesimal generator of the process 
which is defined as : the radial Ornstein-Uhlenbeck process, with dimension 
5 = 2(1 — a), and drift parameter (—//), pushed downwards with parameter 

9 — — . That this infinitesimal generator may be expressed in terms of 
will be discussed after (2.7). 

2.3 We shall now prove a remarkable relationship between the two families 
of diffusions whose infinitesimal generators are found on the RHS of Table 
2. This relationship explains precisely why (Row 1) may be deduced from 
(Row 2), and vice- versa. 

Proposition. 

The following relationship holds with . -6 = — : 

x 2 -«,*ei(t) = x_ h ^ i (4 jT *! Q , Mi H dv) (2.8) 
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meaning that : starting from X = on the LHS, there exists 

(x_„ .^(w), u > J swc/z £/ia£ t/ie relationship (|2.8j) holds. 

Comment about our notation 

In formula (|2.8j) . and possibly several times below, we have written X^ei, 
etc... instead of X i , for some index i. It seemed more appropriate here, 
because of the power 2 on the left-side of (2.7). 

There should be no confusion between the different diffusions X^g and Xi-e^. 
Proof 

We start from the stochastic differential equation satisfied by ^.^(t), 
t > 0) as described (implicitly) in Row 2 of Table 2. 
Then, taking squares, we obtain : 

ft ft K'a / Y 2 \ 

X 2 t =x 2 + 2 / X s d£ s + St + 2 / (/iX s 2 )^ ) ds 

Jo Jo Ka V 2 
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We now define (Y u = Y(u), u > 0) via : 

X 2 = Y^j\ 2 ds^ , (t>0) 



and find that K satisfies : 



r «^ 2+ ^ + (i)ri + u""S(f y -> ds 



(5 (5 1 

since - = — - — , with S = 2 — a = 2(1 — — ), we find that (Y u , u > 0) is 
precisely the diffusion with infinitesimal generator C v }a, with V2u = — , i.e 



We now remark that the proof we have just given for the Proposition re- 
lies upon the identification of the infinitesimal generator of the diffusion 
X_ a ^ (ia)j, as gi ven m Table 2; this identification was obtained from an an- 
alytical identity between Wo t . and K , . (see formula (2.5)). 
We now explain and prove the Proposition without relying on such identities, 
but rather on absolute continuity relationships between the different laws in- 
volved. 
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We now find it a little more convenient to refer to the laws {Q 5 ^} and the 
main absolute continuity result we need is : 

Q}£ = exp {-^{Z t -8t-z)-!j£ Z s d.^ . Q s zWt . (2.9) 

Here, (Z t ,t > 0) denotes the coordinate process on the canonical space 



We now combine this relation 1)2. 9j) with that of the "push downwards" with 
parameter 9, so that, with notations which we shall explain after writing the 
formula : 

Q^ ei (2.10) 



' z \r t 

m(Zi) exp^(9)L t - 9t) exp (~(Z t - 5t - z) - ^ [ dsZ ^j ■ Q 



s 

z \r t 



and we note that for precisely : 9 = — , this relation simplifies as : 

exp (-f z) ipei(z) V 2 

(2.11) 

(The due explanation of the formula (12.101) is that we have combined the 
"push-downwards" formula 1)1. lj) . relative to {Q s z ' M }, - i.e. the function ^ 
and ^(9) are relative to that diffusion - with the preceding formula ()2.9jl ). 

S jj/ 

From now on, we keep : 9 = —. 

2.4 We now consider what becomes of formula [)2.11)1 . once we time change 
both sides with the inverse of I 4 / Z u du, t > I , so that, by a slight abuse 



of notation, the process of reference is now (Z(h), h > 0), with Z defined 
by: 




(2.12) 



Thus, we obtain : 



= t * e W U(9)L u -V-u) .Q% (2.13) 

From the well-known property of time change for Bessel processes (see [12], 
Chap. XI, Prop. 1.11), Q s z is the distribution of a Bessel process of index a/2, 
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i.e. of dimension 5 = 2 — a, that is Q s z = P & z . Again, with obvious notation, 
the right-hand side of (|2.13^ may be written : 



8 * 

and we discover that : 



exp ¥ ^ \L U - ^ . Q 



8 + 

$ (£) = m 



(2.14) 



Again, let us explain, very much in the same spirit, e.g : the first relation : 
(pg_ (z) = e-» g (p 0i {z) in {EH- 

8 

This translates as : 

(exp - ^T (Z)) = eS'E. ( C -^<*>) (2.15) 

where Z simply denotes a BES process with dimension 5, and Z a process 
with law Of'* 4 . This may be well understood by considering the absolute 
continuity relationship (|2.9j) . when we replace t by T (Z). Then, it follows 
from that relationship that : 

Qf'-(exp(-|T (Z))) = e^^/o T °^^ 



<5 .-^r (Z) 



which is precisely (|2.15j) . 

Now, it is well known (see jH], |H, [101) that the Laplace transform of T , 
under the distribution P* of a Bessel process, is given by: 

^ W := £* («p(-£ro)) = 2-*r(f )-' (f )* A', (f ) ■ (2.16) 
Using (I2.14J) . we can recover the expression of (pg± for 9 = — obtained in 



El using the identity El . 
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2.5 We now develop a discussion similar to that made in ()2.8|) . but with the 
downwards arrows j. now changed into upwards arrows f (for the definition 
of these pushed upwards and downwards processes obtained from a diffusion, 
see Pitman- Yor [10J). 
The analogue of formula (|2.1()j) is now : 

Q\T ( 2 - 17 ) 



m{Zt) exp(-et) exp (-|(Z t -5t-z)-!j£ dsz)j . Q 



5 

z \r t 



and we note again that, precisely for : 6 = — , this relation simplifies as : 

- g^fe^ exp (-f /V) . ^ (2.18) 
exp (-f z) ^t(^) \ 2 Jo 1 4 



(We note that this formula is even simpler than 1)2.111) since here there is no 
local time contribution). 

We now continue to develop an analogous discussion to that made in subsec- 
tion (l2~8l) . 

Thus, we time-change both sides of the absolute continuity relation (12.18)1 

with the inverse of 4 / Z u du, t > 0, with Z, as defined from Z in (|2.12)h 

Jo 

We obtain : 



exp ^r(^) / /r >„ 



^ 2 " exp (-§*) ^(*) ^ V § J * Q ^ 
With obvious notation, this right-hand side may be written : 



^ T 0) 

8 1 

with : 

<p e ?Jz) = e-*'<p <n (z) (2.19) 

8 1 

Q z ' ' 2 is the distribution of a Bessel process of dimension 5 = 2 — a with 
drift i.e. 
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The analytical counterpart of (|2.19ft is the companion formula of ()2.6|) (see 
Appendix): 

M ,- f (e) = 4-fr(i-|)Ve/-f (£); 

while the companion formula of (j2. 16j) is: 

• ' (exp -£t,) ) 2 v 2 ; ^ 2 
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a well-known formula which goes back to Kent Pitman- Yor [10J. 

Appendix : On the Whittaker and Bessel - 
Mc Donald functions 

The following formulae involving these classical special functions are found 
in Lebedev 0, to which we refer with numberings such as : (N)* ... 

a) The Whittaker functions M k ^(z) and Wk )fl (z) are a pair of solutions 
of Whittaker's equation : 

„ , ( it, (i-A 2v 



(p. 279,). 

b) Wk,n admits the integral representation : 



r(/i-fc + |)y 

(see Problem 17*, p. 279*). 

c) 

(see Problem 19*, p. 279*). 



_> 



d) In terms of the confluent hypergeometric function there are the 
relations : 

W k ,„(z) = e-f * Q - A; + /i, 2^ + 1; z^j 
12 



(see (9.13.16),, p. 274). 
(see (9.13.15),, p. 274,). 

Taking k = in the above formula for Wk lfl , one recovers c). 

e) In terms of the confluent hypergeometric function <£>, there are the 
relations : 

M k ^(z) = z^ e-4$Q-Jb + /i, 2/i + 1; z^j 
(see (9.13.16),, p. 274,). 

(see (9.13.14),, p. 274,). 

Taking k = in the above formula for M fej/1 , we obtain: 

Mo, M (z) = 4.ro* + i)Vi/,(|) 
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